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Summary 

High-temperature oxygen permeability measurements had determined grain-boundary diffusivities 

(8Dg b ) in bulk poly crystalline alumina in a recent study by Wada, Matsudaira, and Kitaoka. They predict 

-1/ 

that oxygen 8D gb ,o varies with oxygen pressure as P Q /6 and dominates at low ?0 2 whereas aluminum 

+ 3 / 

SZJ gb ,Ai varies with P Q /l6 and dominates at high p 0 2 . In the present study, these relations enable detailed 

evaluations of alumina scale growth in terms of diffusivity and grain size. It is deduced that coupling 
these relations with a modified Wagner treatment for dominant inward oxygen growth produces the 
concise solution: IJ, = k pi -G, = 12 8D gb0 ,mb where IT, is a constant and k pj and G, refer to instantaneous 
values of the scale parabolic growth constant and grain size, respectively. To support this analysis, a 
commercial FeCrAl(Zr) alloy was oxidized at 1100 to 1400 °C to determine k pj , interfacial grain size, IJ, 
and thus SD gbi0 ,mt- The 8D gb , 0 ,,m values predicted from oxidation were 1 0 to 40 times less than those 
obtained by the above solution, but closer than extrapolations from typical high-temperature bulk 
measurements. Furthermore, the activation energy predicted from permeability relations for SD gb 0 ,mt 
(298 kJ/mole) is shown here to be significantly lower than that predicted for bulk alumina (467 kJ/mole) 
under no Po 2 gradient. In comparison, the activation energy determined from oxidation was 375 

kJ/mole. The experimental oxidation results here agree with similar alumina scale studies, especially 
where both k p j and G, were characterized. Temperature-sensitive grain enlargement, equilibrium interface 
P o 2 variation, and grain-boundary diffusivities all affect the Arrhenius behavior of scale growth kinetics. 


1.0 Introduction: Alumina Scale Growth 

Many high-temperature Ni-, Co-, or Fe-based structural alloys and coatings depend on slow-growing 
alumina scales for oxidation resistance. Consequently, alumina scale growth has been the subject of 
intense study over many decades, and many excellent reviews are available that provide both a historical 
perspective of advances and summarize classic key features (Refs. 1 to 4). Briefly, alumina scales are 
shown to grow primarily by oxygen inward grain-boundary diffusion, with some contributions from 
outward aluminum grain-boundary diffusion. The presence of one or more reactive element dopants can 
affect diffusion, most commonly the reduction of outward aluminum diffusion, although no simple 
relationships between diffusivity and dopants have yet emerged. In any event, scale grain size is shown to 
be a critical factor, with most models assuming an inverse dependency between effective diffusion and 
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average grain diameter. Efforts to elucidate diffusional behavior often include l6 0/ lx 0 two-stage or 
double oxidation tracer studies with nuclear reaction analysis (NRA) or secondary ion mass spectroscopy 
(SIMS) reconstruction of diffusion profiles. Others utilize growth rate or electrical conductivity variations 
with P o 2 , then model the ionic transport contributions as a function of the Po 2 gradient (Refs. 5 to 1 1). 

This latter group provided the relations necessary for describing scale growth according to the Wagner 
relation in which diffusivity is characterized as a function of the chemical potential gradient across the scale. 
Others address the Po 2 gradient in the integration limits for the modified Wagner equation in which Ain 

Po 2 becomes the salient factor (Refs. 2 and 12 to 15), and many of these will be highlighted later. 

The present study is not a critical review or synthesis of the various kinetic studies. Rather, it applies 
grain-boundary diffusivity relations recently obtained via permeability measurements to the case of scale 
growth. In those studies diffusion in bulk alumina samples had been elegantly determined through oxygen 
permeability measurements as a function of Po 2 gradients (Refs. 16 to 19). An oxygen flux had been 

produced across a thin polycrystal of alumina by establishing a Po 2 gradient using gas tight seals and 
controlled 0.01 to 10 vol% CT/Ar or 0.01 to 1.0 vol% tE/Ar atmospheres. Typically, a high- P Qi regime 

was studied by keeping one side of the membrane at 1 Pa O 2 and varying the other to 1 0 5 Pa. 

Alternatively, a low- P 0 , regime was explored by keeping one side of the membrane at 1 Pa O 2 and 

varying the other down to 10’ 9 Pa. In some cases the entire pressure regime was studied as a closer 
simulation to a growing scale. Experiments were performed over the 1500 to 1700 °C temperature range, 
and O 2 pressure was monitored via zirconia electrochemical cells. Oxygen permeability was determined 
from the change in steady-state P Qi on one side as the P Qi was changed on the other. The measured 

permeability and flux were converted to oxygen and aluminum diffusivities at the low- and high- p 0 2 

-1/ +3/ 

regimes, respectively, as they followed the expected P. ' /6 and R/ 16 dependencies for n-type oxygen 

vacancy Vo and p-type aluminum vacancy V A i defect equilibria. 

Experimentally, the present study draws from early grain size measurements made for alumina scales 
formed on a commercial FeCrAl(Zr) alloy, Hoskins 875 (Hoskins Manufacturing Co.), originally 
intended to verify grain size effects on oxidation rates in early transmission electron microscopy (TEM) 
studies of NiCrAl scales (Refs. 20 and 21). However, the primary purpose at present is to reconcile 
alumina growth kinetics and grain-boundary diffusivity obtained by combining the permeability equations 
with the Wagner model. To accomplish this, a number of relations must be derived. In Section 2.0 we 
present the basic grain-boundary diffusion equations resulting from the permeability studies. Their 
incorporation into the Wagner integral, using the concept of effective diffusivity D eS and grain size, yields 
the general interrelation of the parabolic growth rate constant k p and oxygen grain-boundary diffusivity at 
the interface §/J g b.o.im- In Section 3.0, the experimental data set of k p and grain size are compared with 
values predicted from 5D gb- o. inl . Since the latter is dependent on the p o 2 at the metal-scale interface, the 

thermodynamic calculation of P 0 , e q > s first presented. The concept of a constant oxidation product, n, is 

then introduced and used to extract SZ? g b,o,mt from experimental values of k p and grain size. Comparable 
values and activation energies from the literature are noted throughout. Finally, grain-boundary 
diffusivities obtained from permeability, oxidation, and ls O tracer studies of scales and bulk alumina are 
compared and discussed. 

A list of symbols used in this report is given in Appendix A to aid the reader. 
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2.0 Development of Diffusivity and Oxidation Equations 

This section addresses new relations for A1 and O grain boundary diffusivity as obtained from oxygen 
permeability studies. These are then combined with the standard Wagner analyses that relates scale 
growth rates to diffusivities. 


2.1 Permeability and Grain-Boundary Diffusivity 


The assumption of the primary ionic species is essential to the treatment of diffusion across a P o 2 

and [V&], 


Vo 


gradient in alumina. For the sake of brevity, we consider charged O and A1 vacancies ( [ 
respectively) as likely candidates to be formed at interfaces and incorporated as charged polarons 
diffusing within grain boundaries (Ref. 22). The standard relations connecting defect concentrations with 
Pq are developed in Appendix B. These relate the reaction of A1 and O to form A1 2 0 3 and the formation 


of charged 


Vo 


and [FI)] in accordance with equilibrium thermodynamics using Kroger-Vink notation 
(Ref. 23). Then, electroneutrality is defined for the primary charged ionic vacancies and electron or hole 
charge carriers. By combining the latter with the equilibrium constants and solving for defect 

t~Y' 

' 0 2 


concentrations, one obtains 


Vo 


proportional to P. ^ at the low- P ()n surface (n-type) and \V t 


All 


proportional to P n 


at the high- P 0l surface (p-type). Indeed, the grain-boundary diffusivities calculated 


from permeability studies accordingly yield expected P () '' (Eq. (1)) and P () 16 (Eq. (2)) dependencies at 
low (<1 Pa) and high (>1 Pa) p 0 2 , for oxygen and aluminum vacancies, respectively (Ref. 24). Fe, Ni, Y, 

Elf, and Zr dopant effects as well as Schottky and Frenkel defect equilibria are significant considerations 
and have been considered in numerous permeability, conductivity, and oxidation studies. Flowever, the 

-1/ + 3 / 

present discussion will refer to the simple P Q /6 and P () u> relations in Appendix B, since they agree with 

the permeablility results for grain-boundary diffusion in pure alumina. Resulting from this are highly 
useful relations for both O and A1 grain-boundary diffusion (Eqs. (1) and (2), respectively) as a function 
of temperature and oxygen pressure. From Reference 24, 


5T) gb o = 2.207 x 10 9 exp 


8D 


gb,Al 


: 2.475 xl0“ 5 exp 


-467UVK 

m 3 

(1) 

RT 

J 0 2 ,I 

s • Pa Vi 

-604kjV +X( 

m 3 

(2) 

RT 

/o 2 ,n 

s-Pa + ^ 6 


where R is the ideal gas constant, T is temperature, I refers to the low- p o 2 surface, and II refers to the 

high- P Qi surface. These relations therefore define 8D gb 0 and 5D gh . A | as a function of temperature and 

oxygen pressure, allowing full predictions of transport across a Po 2 gradient, such as in a growing 

alumina scale. These authors (Ref. 24) applied such equations to FeCrAl oxidation literature values and 
provided a basis for the extended comparisons developed in this paper. Combined with the appropriate 
chemical potential gradients, these diffusivities present a strong case in support of primary inward oxygen 
flux and secondary aluminum outward flux for alumina scale growth (Ref. 24). 

As the following sections rely on a number of derivations from various kinetic relations, it may be 
useful to refer to the flow diagram in Figure 1 for an overall perspective of corresponding sections. 
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Equation flow diagram 



k p = parabolic growth constant 

SDgb = grain boundary diffusion product 

k p j = instantaneous parabolic 

growth constant 

Peq = equilibrium pressure 

k p o = parabolic growth constant of oxygen 

D e ff o = effective oxygen diffusivity 

Dgb.o = g ra i n boundary diffusivity for oxygen 


Calculate k p o (from D e ff o ~ 25Dgb,o/G) 
in Wagner integral; Secs. 2.2, 3.5 


Calculate 5 D g b (O, Al versus T, P eq ) 
permeability relations; Sec. 3.1 



SOgb.o.int 

Gi 

W-W-S 


grain boundary diffusion product of 
oxygen at scale-metal interface 

instantaneous grain size 
Wagner-Wada-Smialek 


Figure 1. — Derivation scheme relating oxidation to diffusion, then comparative analyses of predicted and measured 
oxygen diffusivities. 


2.2 Oxidation Models 

Oxidation processes are generally described by the Wagner relation for the counterdiffusion of ionic 
species across the chemical potential gradient in the scale, with z a and z c ionic charge and D„ and D, 
diffusivity of each anion (a) or cation (c) species (Ref. 23). With the parabolic rate constant k p defined as 
given below in (Eq. (3)), the appropriate form of the Wagner relation follows as (Eq. (4)). (Note: A factor of 
1/2 is often employed before the integral, but corresponds to the less frequently used differential form of k p 
= lx-cbc/df, rather than the k p = 2x-dx/dt definition used here; note t is time. Also, for the purpose of direct 
comparisons to diffusivities, the thickness form of scaling constants, in m 2 /s, was ultimately adopted.) 


x 2 = kpt ; 

k v j-2x — 
P ’ 1 d t 

(3) 

A), .gas r 

) 


J 7 

D c +D a >dln Pq 2 

(4) 

P l z a 

mD 2 ,int 

J 



where z is the valence of anions and cations in the compound. Now, turning to fine grain alumina scales, 
the importance of grain-boundary diffusion is well recognized and illustrated by the Elart/Smeltzer 
equations for alumina in Equations (5) and (6) (Refs. 25 to 27). Elere the effective diffusion coefficient 
D efi is a function of the area fraction of short-circuit grain-boundary paths/ often simplified as 2S/G, where 
8 is the grain-boundary width and G is the grain size (diameter). In the case of alumina, it is generally 
observed that D sh » D L (lattice diffusivity), leading to the following simplification in Equation (6): 

D cff =(l-/)D L+ /D gb (5) 
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D eff ~ fD gb 


( 6 ) 


25P gb 
G 

Most transport marker and tracer studies have indicated that the primary growth of alumina scales 
takes place by oxygen inward diffusion, although some outward growth of alumina scales has been 
clearly documented (Refs. 28 to 3 1). For the purpose of illustration, it is assumed here that D g b ,o » F) gbjA i; 
thus the latter will be omitted from the Wagner integral at present. Flowever, it is also well established 
that grain “growth” can occur, increasing with scale thickness and temperature. Thus the scale growth rate 
becomes grain-size dependent and is best treated as an instantaneous value (subscript The resulting 
expression for k pJ is obtained by substitutions in Equations (3) to (6) to yield 


k Pti = 2x— = 


dv 

At 


0,,gas 


- I 


25 D 


gb.O 


din R 


O, 


(7) 


If the diffusion product 8D g b,o and grain size are arbitrarily assumed to be invariant across the P 0 
gradient, the following oft-used integration for k px results: 


k 


PA 


2dPgb,0 

Gi 


Mn Pq 2 


[LNP] 


( 8 ) 


Flere 8 Z)* b 0 is designated as an approximate or average value and is referred to as the “LNP” solution. 

However, with the complete Pq 2 dependence provided by the permeability studies of Wada, Matsudaira, 

and Kitaoka (Ref. 24), the more precise relations for k pJ are obtained by integrating below, where A is the 
pre-exponential constant in Equation (10): 


k 


pj 


db.gas 2d exp 

I — 


G, 


r O, ,int 




d P, 


o-, 


(9) 


k P,i = 


>(2.207 x 1 O' 9 )exp(- 467 W 


Gt 


rt, 6a p-y* 
u 2 


128 D 


gb.O.int 

~Gi 


[W-W-S] 


( 10 ) 


a p 76 ~ _p 76 

Finally, since o 2 ~ 0 2 ,int , it is seen by inspection that the Wagner integral of the Wada, 

Matsudaira, and Kitaoka diffusivity k pJ is reduced to a quite simple relation in 5D gb 0 , abbreviated here as 
the Wagner- Wada- Smialek (W-W-S) solution (Eq. (1)). Furthermore, the deviation in 5/J gb ,o calculated 
from the W-W-S solution can be shown as an offset from the constant 8 0 (LNP solution, Eq. (8)) by 
the factor A In Pq 2 /6 in the final relation: 


k P J = ' 


128 D, 


gb,0 


Gi 


7 5f>gb ; o A In P Ql 
G, 


(11) 
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ZA 111 I f) 

5Z) gb>0 [W-W-S] - 6 ° 2 SP g \ Q [LNP] (12) 

3.0 Experimental and Predicted k p and 5Z) gb ,o 

Given these derivations and expressions, predictions for k p (W-W-S) from the Wagner integral will be 
compared, using these permeability expressions for grain-boundary diffusivity, with those determined 
experimentally. The results for the reverse calculation, showing what values of 5/J„b.o are obtained from 
applying the experimental oxidation results to the W-W-S solution, are then compared to those defined by 

the permeability expression. Equation (1). Results will also be presented for SD gb ,o > reflecting the 
departure that occurs if one assumes diffusivity is a constant or average across the P 0 , gradient in the 
scale. This begins with the necessary thermodynamic treatment to obtain P 0i int for the given alloy 

versus temperature, allowing the calculation of 8D gh 0 from the permeability expressions. This is followed 
by a summary of the experimental oxidation results needed to make these comparisons. The k p and 8D gb 0 
determined from experimental oxidation and diffusion models can then be compared to those predicted by 
the permeability relations. 


3.1 Calculations Versus Temperature 

The Wagner integration is performed over the P Qi limits for the scale, from the equilibrium 

dissociation pressure at the oxide-metal interface to that of the external gas. The former is obtained from 
the equilibrium constant of the oxidation reaction and the activity of aluminum in the alloy. The approach 
follows previous studies that predict P Qi int eq from thermodynamic databases and estimate diffusivity 

from scale growth kinetics (Refs. 10, 12, 14, 15, and 32). For the present example, using a typical 
Hoskins 875 heater alloy, approximately Fe-22.5Cr-5.5Al-0.5Si-0.2Zr by weight, A1 activities were 
calculated by using the Pandat thermodynamic software package (CompuTherm LLC, Ref. 33) in 
conjunction with the Paniron database and equilibrium constants obtained from the FactSage 
thermodynamic database for chemical reactions (Center for Research in Computational Thermochemistry, 
Ref. 34). The pertinent relations are given below as Equations (13) and (14). The obtained activities a A i, 
equilibrium constants K eq , and P ()n eq are given as a function of temperature in Table 1 and Figure 2. 

Given the high degree of linearity, a relation for the equilibrium pressure is also presented (Eq. (15)), 
obtained by regression of the data with a correlation coefficient of 1.000. This enabled the corresponding 

A P n /6 and A 1 n P () to be evaluated, as listed in the table. (These essentially converge to - P r . /6 . and 

-In P Q int , respectively, because of the extremely low values of P Qi int compared to P Qi gas .) 


2A1 + X C) 2 = A1 2 0 3 

(13) 

K - 1 

iv eq 9 3/ 

(14) 

£ i n i ni5 f-1012kJ^| 

= 6.117x10 exp Pa 

l RT J 

(15) 


NASA/TM— 2013-217855 


6 



TABLE I.— TEMPERATURE EFFECTS ON THERMODYNAMIC FACTORS FOR ALUMINA SCALES IN 


EQUILIBRIUM WITH FeCrAl(Zr) a b 
[Oxygen pressure at gas surface -Po 2 ,gas = 2.027 x 10 4 Pa.] 


Temperature, 

T, 

°C 

Aluminum 

activity, 

Oai 

Equilibrium 

constant, 3 

K qq 

Oxygen pressure at 
metal interface, 

p 

1 02,int 9 

Pa 

, 

Pa~^ 

AlnP 0 , , 
Pa 

1000 

9.760xl0“ 04 

1.400xl0 52 

1.778xl0" 26 

19573.36 

69.21 

1100 

1.747xl0“° 3 

1.245xl0 47 

1.906xl0" 23 

6118.07 

62.23 

1200 

2.896xlO“ 03 

5.409xl0 42 

7.861xl0" 21 

2242.40 

56.21 

1300 

4.509x1 0“° 3 

8.490xl0 38 

1.497xl0“ 18 

934.75 

50.96 

1400 

6.665xl0“ 03 

3.822xl0 35 

1.514xl0" 16 

432.97 

46.34 


K cq - i/[ahp £ ) for 2A1 + y 2 o 2 = Al2 o 3 


b Wagner integral limits on P 02 e!| fixed by corresponding interfacial aluminum activity and equilibrium constant. 



Figure 2. — Equilibrium interfacial Pq 2 jnt calculated from thermodynamics 
of FeCrAI(Zr)-scale interface using aluminum activity and its reaction with 
oxygen. Slope of log P versus 1/7" found to be equivalent to -1012 kJ/mole. 


At this point, the relationships between k p , 8D„ bi0 and ^o 2 ,eq can now be employed to evaluate 

experimental data for scale growth and grain size as a function of time and temperature. However, it is 
instructive to first illustrate some broad implications for 8/X b-0 and 8D gb , A i at both sides of the scale and 
for bulk alumina as well. The gas surface values are obtained from Equations (1) and (2), using ^) 2 ,gas = 

2.027xl0 4 Pa (0.2 atm). However, the interfacial equilibrium ^o 2) mt (i-e-, referring to the metal-oxide 
interface) is given by Equation (15). Both result in a modified A\ exp (Q'JRT ) form, where i designates 
an O or A1 ion. A- gas is given by A i P 0 2 , gas and Q\ = Q, from Equations (1) and (2). A\ int is given by 
A i AP o 2 ,gas and 2/, int -Q + n iQp 0 , i as summarized in Table II (and compared further in Section 5.0). 
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TABLE II.— GRAIN-BOUNDARY DIFFUSION COEFFICIENTS 5 Z) gb , x PREDICTED FROM 

PERMEABILITY RELATIONS 3 


Temperature, 

T, 

°C 

Grain-boundary diffusivity, 
m 3 /s 

Oxygen 

Aluminum 

Gas surface 
(bulk), 

^*gb,gas,0 

Metal 

interface, 

^gb.int.O 

Gas surface 
(bulk), 

^^gb.gasAl 

Metal 

interface, 

80 g hj n,^A, 

1000 

2.912xl0“ 29 

2.976xl0“ 24 

2.618xl0“ 29 

6.057xl0“ 35 

1100 

7.241 xlO -28 

2.313xl0“ 23 

1.671xl0“ 27 

1.431xl0“ 32 

1200 

1.164xl0“ 26 

1.363xl0“ 22 

6.066x10 _2S 

1.606xl0“ 30 

1300 

1.314xl0“ 25 

6.416xl0“ 22 

1.395xl0“ 24 

9.883xl0“ 29 

1400 

l.lllxl 0 -24 

2.512xl0“ 21 

2.205xl0 -23 

3.712xl0“ 27 


Pressure exponent, n 

-1/6 

-1/6 

+3/16 

+3/16 

Pre-exponential factor, A\ m 3 /s 

4.23xlO~ 10 

5.16xl0“ 12 

1.59xl0"° 4 

2.26xl0“° 2 

Activation energy, 3 Q\ kj/mole 

467 

298 

604 

794 


a Pre-exponential factors \A'.j and activation energies \Q'.j are adjusted by pressure, exponents, and energy associated with 


P 02 ,eq in the general Arrhenius form of 5 D gb l = A- exp(— i Q- /RT ) . 



Oxidation temperature, 1/T, 10 4 K 1 

Figure 3. — Grain-boundary diffusivities predicted for alumina scales on 
FeCrAI(Zr) using permeability relations for aluminum and oxygen at the 
gas surface and scale-metal interface. Magnitude and slopes change 
because of large Po„ gradient in scales and temperature-sensitive 

0 3/16 n -1/6 , f 

Po 2 ,int or Po 2 ,int factors. 


The limiting relations for 8 D„h at each reaction surface for each ion are plotted in Figure 3. The 
highest 8/Jgb is predicted for oxygen at the scale-metal interface, consistent with primary inward growth 
of alumina scales, as confirmed by virtually every diffusion study for a-AFC^ scales. Oxygen diffusivity 
at the interface is about 4 to 5 orders of magnitude higher than at the gas surface. Conversely, the highest 
SDgb for A1 occurs at the outer gas surface. It is 3 to 5 orders of magnitude higher than at the scale-metal 
interface and similar to or slightly above the corresponding 8T> g b,o,gas- 
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The broad arrows in Figure 3 indicate the changes in diffusivity obtained as the scale is traversed 
inward. The transition between predominantly A1 outward to predominantly O inward has been predicted 
by Wada, Matsudaira, and Kitaoka (Ref. 24) and Matsudaira et al. (Ref. 18) to be extremely close to the 
gas surface, at just 0.01 percent of the thickness, consistent with the frequent observation of some slight 
outward growth of alumina scales. Thus they project that double oxidation experiment utilizing l6 0/ 18 0 
diffusion profiles actually ’’probe” levels deeper than this transition and reflect transport still dominated 
by 8£> gb ,o. 

Finally, it should be pointed out that 8/J gh . 0 measured on bulk alumina with no P 0i gradient (i.e., 

self-diffusion) should follow the relations for SD g b,o, g as, as appropriate to the external pressure. This is 
suggested as “bulk” behavior in Figure 3 and would be quite lower than (chemical) diffusivity throughout 
the scale (arrows). 


3.2 Oxidation Procedures 

To provide comparisons between permeability predictions for scales and actual data, we draw from an 
older dataset, originally intended to augment some early scale diffusivity discussions for NiCrAl 
(Refs. 20 and 21). For this purpose, a commercial ferritic FeCrAl(Zr) heater alloy, Ffoskin s 875, was 
utilized. The composition is listed in Table 111. The alloy was obtained as rolled strip ~1 mm (~40 mil) 
thick and sectioned into 1.25- by 2.54-cm (0.5- by 1-in.) coupons. These were metallographically 
polished on emery paper down to a 600-grit finish, washed in detergent, and rinsed ultrasonically in 
ethanol. Samples were oxidized in high purity alumina boats at 1100, 1200, 1300, and 1400 °C in static 
ambient air. Weight change was measured for samples pulled at various times up to 1000 h (at 1100 and 
1200 °C), 500 h (at 1300 °C), and 100 h (at 1400 °C). (Longer tests at the higher temperatures resulted in 
chemical failure of the scale at temperature.) Thus, each weight for each time and temperature exposure 
corresponds to a specific sample. As a consequence, these results may be subject to sample-to-sample 
variations, especially regarding dx/dt. In a few cases minor spalling was measured in the small crucibles 
containing the samples and was incorporated into the total mass gain, but it made little difference. 

Samples were pulled for grain-size measurements after oxidation for 0.1, 1, 10, 50, 100, 200, 500, and 
1000 h, where applicable. Grain size was measured at the oxide-metal interface by detaching the scale. To 
accomplish this, aluminum stubs were attached to the scale with epoxy, then detached by pulling in a 
tensile adhesion test jig. Numerous representative micrographs of the underside of the alumina grains 
were obtained by SEM. Also, the corresponding grain imprints in the exposed metal side of the coupon 
were measured. The recorded grain size was the largest diameter of a grain, averaged over -100 grains 
measured for a given test condition. 

The underside of representative scale samples are presented in Figure 4. Figure 4(a) shows the 
noticeable temperature effect on the enlargement of grains formed after 100 h, and Figure 4(b) shows the 
time effect for grains formed at 1400 °C. At higher temperatures, a tendency for dimpled grains was 
noticed. These were matched and in contact with intergranular protrusions on the metal side and so were 
not voids. These features may indicate an increasing inward oxygen flux in the vicinity of the grain- 
boundary relative to that along the lateral interface at higher temperatures. 


TABLE III.— COMPOSITION (wt.%) OF FeCrAl(X) ALLOYS WHERE OXYGEN GRAIN-BOUNDARY 
DIFFUSIVITY 5D gb 0 IS PREDICTED FROM OXIDATION KINETICS AND GRAIN SIZE 


Study 

Alloy 

Alloy composition, 
wt% 

Fe 

Cr 

Al 

Si 

Zr 

Y 

Other 

Smialek (current study) 

Hoskins 875 

Bal. 

22.5 

5.5 

0.5 

0.2 

.... 

.... 

Naumenko (Ref. 35) 

FeCrAlY 

Bal. 

19.9 

5.00 

0 

.... 

0.051 

.... 

Messaoudi (Ref. 36) 

Imphy 

Bal. 

22.6 

4.36 

0.85 

0.13 

.... 

0.38 Mn, 0.35 Ni 

Ramanarayanan (Ref. 10) 

MA956 

Bal. 

20.0 

4.5 

.... 

.... 

.... 

0.5 Ti, 0.5 Y 2 0 3 
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G- 0.7 pm 1.0 |jm 2.2 |jm 5.5 pm 



t= 1 h 1.0 h lOh lOOh 

(b) 


Figure 4. — Representative scale (underside) structures stripped from oxidized FeCrAI(Zr) alloy, showing alumina grain size G 
increases from 0.2 to 6 pm (on average), (a) Oxidation at 100 h for different temperatures T. (b) Oxidation at 1400 °C for 
different times t. 


3.3 Experimental Weight Gain and Grain Size Data 

The weight gain and interfacial grain size data are given in Table IV. Although the details of the 
kinetic behavior are the subject of another study, it is noted that the weight gain was subparabolic, 
yielding a time exponent oft , while grain size followed a ~t ' dependence, on average. This is quite 
consistent with D s{f being inversely related to grain size, where a f 0 375 oxidation law was predicted for the 
t 0 25 grain “growth.” (Oxidation time exponents (\/n) and grain growth time exponents ( 1/m) are related 
by n = 2m/(m— 1), Refs. 21, 37, and 38). The t° 37 rate is similar to the cubic FeCrAl(Y) oxidation law 
presented by Quadakkers et al. (Refs. 38 and 39). 

Similary, the grain expansion can be assessed in relation to the findings of Naumenko et al. for the 
case of FeCrAlY oxidized at 1200 and 1250 °C (Ref. 35). Here grain size G, elegantly measured by 
EBSD (electron backscatter diffraction) of FIB (focused ion beam) cross sections, was found to be a 
temperature-invariant function of scale thickness x in the form G = 0.27 + 0.1 3x (both in pm). Using 
weight gain obtained in the present study of FeCrAl(Zr) (converted to 5.34 pm thickness for each mg/cm 2 
of oxygen gained) from Table IV, interfacial grain size versus scale thickness was plotted in Figure 5. It is 
seen that the thinner scales cluster along the single curve of the Naumenko relation; however, positive 
deviation appears significant at the higher 1300 and 1400 °C temperatures. In any event, the importance 
of increasing grain size for thicker scales is clearly demonstrated. This treatment simply employs the 
instantaneous grain size (Table IV). 


NASA/TM— 2013-217855 


10 



TABLE IV.— PERTINENT SCALE GROWTH AND GRAIN SIZE DATA 
USED TO ESTIMATE IT, AND 5D gb>0 FOR HOSKINS 875 a 


(a) 1100 °C 


Time, 

h 

Scale weigh gain, 

x, 

mg/cm 2 

Scale growth rate, 
dx/dt 
mg/cm 2 -h 

Grain size, 
G, 
pm 

Oxidation product, 3 

n, 

pm 3 /h 

0.1 





1 

0.027 


0.56 


10 

0.192 

1.02xl0“° 2 

0.42 

4.69xlO“ 02 

50 

0.438 

4.25xlO“ 03 

0.59 

6.28xlO“ 02 

100 

0.591 

2.24x1 0 -03 

0.70 

5.28xlO' 02 

200 

0.752 

1.36xlO~ 03 

0.62 

3.61xlO' 02 

500 

1.010 

3.92xlO“ 04 

1.01 

2.28xlO“ 02 

1000 

1.077 

4.80x1 0“ 05 

1.18 

3.48xlO“ 03 


(b) 1200 °C 


0.1 





1 

0.152 




10 

0.453 

1.66xlO“ 02 

0.53 

2.27xlO -01 

50 

0.731 

6.67xl0~ 03 

0.80 

2.22xlO“ 01 

100 

1.153 

3.94xlO“ 03 

1.05 

2.72xlO“ 01 

200 

1.38 

2.10xl0~° 3 

1.07 

1.77xlO“ 01 

500 

1.725 

8.00xl0~° 4 

2.18 

1.72xlO“ 01 

1000 

1.979 

4.19xlO~ 04 

1.66 

7.85xlO' 02 


(c) 1300 °C- 


0.1 

0.125 

9.50xl0 -01 

0.45 

3.05 

1 

0.492 

1.87xlO~ 01 

0.75 

3.93 

10 

0.910 

3.62xlO“ 02 

1.25 

2.35 

50 





100 

2.147 

8.50xl0~° 3 

2.22 

2.31 

200 

2.801 

5.17xlO“ 03 

3.67 

3.03 

500 

3.713 

1.20xl0~° 3 

4.09 

1.04 

1000 






(d) 1400 °C 


0.1 

0.246 

1.70xl0 +o ° 

0.58 

1.38X10 401 

1 

0.764 

3.43xlO“ 01 

1.27 

1.90X10 401 

10 

1.731 

7.12xlO~ 02 

2.31 

1.62X10 401 

50 





100 

4.063 

1.13xlO“ 02 

5.51 

1.44X10 401 

200 





500 





1000 






Hi, is oxidation product n, (= k p /G h where k pi is the instantaneous value of the 
parabolic growth constant and G, the instantaneous grain size) at the 
low-oxygen-pressure surface. 5D gb 0 is the oxygen grain-boundary diffusion product. 


NASA/TM— 2013-217855 


11 







Figure 5. — Correlation of interfacial scale grain size G with scale 
thickness x estimated from oxidative mass gain from interrupted 
FeCrAI(Zr) oxidation at 1 100 to 1400 °C. Comparison to relation 
for scale cross sections produced by Naumenko et al. for 1200 
and 1250 °C oxidation (Ref. 35). Here C is 0.2273 pm 3 /h and 
t is time (5.34 pm thickness conversion for every 1 mg/cm 2 
oxidative mass gain). 



Figure 6. — Relationship of FeCrAI(Zr) instantaneous parabolic 
oxidation constant kpj\o grain size G. Initial portions (at ~0.5 pm) 
indicate an inverse G" 1 relationship holds at all four temperatures 
(1 mg 2 /cm 4 *h = 7.918x10~ 15 m 2 /s). 

The corresponding specific, individually measured, instantaneous k p j were determined from 
Equation (3), obtaining dx/dt graphically from the slope of x versus t plots. The variation with grain size is 
given in Figure 6. The initial values of k p i decrease at the predicted rate according to G '. These portions 
are thus consistent with the treatment of k pj as a function of D e ff=28D gb , 0 /Gi (Eq. (10)), with 8Dg b . 0 
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assumed to be constant with time. They are therefore considered suitable for assessing grain-boundary 
diffusion as an invariant parameter in the scaling process. The minor fluctuations in each curve are 
believed to arise from graphically determining 2x-dx/df from a single oxidation curve generated by 
multiple samples, with perhaps some sample-to-sample variation, rather than using one continuous curve 
for one sample. Additionally, the variation in grain size could be substantial within any given sample, but 
the average values of many measurements should provide valid trends. 

3.4 Arrhenius Behavior of k pi 

For purposes of consistent comparisons and subsequent calculation of a deconvoluted activation 
energy, the instantaneous k pp were determined at the same grain size (0.5 pm) for each temperature. Other 
grain sizes would also be appropriate for this analysis, up to about 1 pm, according to Figure 6. An 
Arrhenius plot of the aforementioned k p ^ 0.5, is presented in Figure 7. The fitted relation for k p ^ 0.5, 

(r 2 correlation coefficient of 0.98) is given by Equation (16) and shows an activation energy of 
375 kJ/mole. 


' : p,/,cxptl 


. „_ 3 [ -375kJ 

= 4.29x10 exp 

1 RT 


m 2 /s 


(16) 


This curve is now compared with data from the literature, where grain size was provided in addition 
to kinetics. Both parameters are available unequivocally from the study of Naumenko et al. (Ref. 35). 
Flere, well-defined data were examined to determine the equivalent thickness for a 0.5-pm grain size for 
scales grown on FeCrAlY (Table III) at 1200 °C. Their equation for grain size yields a scale thickness of 
1.77 pm for an interface grain size of 0.5 pm. This was then used in their cubic polynomial describing 
scale growth, 0.13x 3 /3 + 0.27x 2 /2 = Ct, where C can be shown to be 0.2273 pmVh. Differentiating and 
solving for 2x dx/df yields 2 C/G. At G = 0.5 pm, 2 C/G yields 0.909 pm 2 /h or 2.525xl0~ 16 m 2 /s for k p j. 
The latter value, plotted as the square symbol in Figure 7, is seen to fall directly on the 0.5-pm regression 
curve for the present FeCrAl(Zr) alloy. This agreement provides some degree of confidence with the 
present technique in which separate samples were used to determine k p j and G,. 



Oxidation temperature, 1/T, 10 4 K -1 

Figure 7. — Instantaneous k p j, measured for G, = 0.5 pm, compared to 
predictions from permeability relations. Experimental data at G = 0.5 pm 
all agree for FeCrAI alloys studied. 
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Similarly, the oxidation behavior of an Imphy FeCrAl(Zr,Si) alloy (Imphy Alloys, Inc.) 

(see Table 111) by Messaoudi, Huntz, and Lesage (Ref. 36) is presented as open diamonds. Here the 48-h 
average parabolic growth rate kp Avg was reported, but it corresponds to similar but slightly varying G 48h = 
0.3, 0.6, and 0.8 pm at 1000, 1 100, and 1200 °C, respectively. Again, the data are seen to fall near the 
present data or regression line. Finally, the average 100-h data of Ramanarayanan et al. (Ram) from 
Reference 10 for MA956 (an FeCrAl-Y 2 03 alloy, Special Metals Corp.), Table III, oxidized at 900 to 
1200 °C, are presented according their Arrhenius relation, converted to kilojoules and meters squared per 
second: 


k p i ioo [Ram] = 1 .20 x 1 0 2 exp 


-388 kJ 
RT 


2 

m /s 


(17) 


Here it should be noted that a relatively constant grain size of 0.5 pm was reported for the 100-h TGA 
(thermogravimetric analysis) test (increasing to only 0.75 and 1.25 pm for 300- and 1000-h exposures, 
respectively). This effectively allows direct comparison with the present k pJ data and is seen to lie 
essentially on the same line. 


3.5 k p j Predicted From W-W-S Permeability Relations 

Having established a valid Arrhenius plot of k p j (i.e., without grain growth), it is now possible to 
compare predicted k p from the permeability studies and the W-W-S solution (Eq. (10)) developed above, 
all for a grain size of 0.5 pm. Using the permeability equation reported for 6D g b,o and the values for 

A P n /( ' from Table I, the values shown as “x” in Figure 7 are produced. These are higher than those 

determined from oxidation experiments and result in an activation energy of 298 kJ/mole, somewhat 
below the 375 kJ/mole determined from oxidation. The Arrhenius dependency obtained with r 2 = 1.000, 
when combining equations for dD gb0 and P 0n int in the W-W-S solution (Eqs. (1), (10), and (15)), is 

given by 


k p i [W-W-S] - 1.24 xlO -4 exp 


298 kJ 
RT 


m 2 /s 


(18) 


As will be discussed later, the discrepancy between the refined oxidation data (Eq. (16)) and the derived 
W-W-S solution (Eq. (18)) may be due in large part to the fact that the former pertains to doped alumina 
scales while the latter pertains to pure alumina. 


3.6 k p j Prediction Assuming P Qi -Independent Diffusivity (LNP Integral) 

As will be discussed below, analyses of alumina scale growth frequently assume a constant diffusivity 
across the scale. The permeability studies of Matsudaira (Ref. 18) suggest that this is somewhat true for 
much of the thickness for the proposed oxidation case modeled for FeCrAlY (Y2O3) alloys at 1100 °C and 
an interfacial P 0i of ~10 23 Pa. Their model predicts 8D gb across the much of scale is close to that given 

by the 7 q 2 at the interface (10 23 Pa). A simple comparison of the LNP solution to the W-W-S solution 
was given as Equation (12) and can be shown to be approximately a factor of 8 to 11 higher: 


k p [LNP ] _ AlnP 02 lnP 02 int 

k p [W-W-S] 6 ~ 6 


(12b) 
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This is presented as the topmost LNP curve in Figure 7 and is seen to have a slightly lower activation 
energy and to deviate further from the experimental data. Alternatively, experimental values of k p>i can be 
used to predict the grain-boundary diffusion product 5Z) gb ,o,mt by both the W-W-S or LNP approaches, 
then compared to those predicted by permeability studies from Equation (1) as well as values obtained 
from tracer studies on bulk alumina. These results are presented in the following section. 

4.0 8Z)g b ,o Predicted From Permeability Relations and Oxidation Data 

The permeability studies have offered relations allowing for the prediction of 8D gb ,o once the P 0 int 

has been determined. This enabled the prediction of oxidation rates for given scale grain sizes using both the 
LNP and W-W-S solutions to the Wagner equation (Eqs. (8) and (10)) and comparison to those actually 
measured. In the following the reverse occurs; namely, estimates of SD^, o are obtained from the oxidation 
data and then compared to those predicted by the permeability and 1S 0 tracer studies of bulk alumina. 

4.1 Time-Invariant Oxidation Product, II 

In order to properly assess the calculation of 8/Ab.o.mt from oxidation data, it is first useful to 
introduce the concept of a time-invariant constant scaling product, IT, = k p j- G,. This follows from 
Equation (10) when the inputs to the Wagner integral, the diffusion product SZ) gb and P 0l at the interface 

(fixed by interface composition), do not change with time. Accordingly, the experimental values of IT, are 
presented as a function of time for the four temperatures of this study in Figure 8. (The use of t m as the 
x-axis is of no mechanistic significance; it was adopted merely to provide a reasonable comparison of the 
data without going to the extreme of a log scale). The plot does suggest invariance of IT,, at least for most 
of the initial portion of the oxidation curves. The dashed lines indicate average values for the invariant 
regions (first three or four data points to about 200 h). An appreciable drop in II, is noted at long 500- to 
1000-h exposures and suggests a reduction in 8D gb or an increase in / J (); int , perhaps due to aluminum 

depletion. 



Figure 8. — Invariance of instantaneous oxidation product, n,, for 
FeCrAI(Zr) reflecting growth by constant 6D g b,o- Comparison to 
invariant 1200 °C values determined for FeCrAlY by Naumenko 
et al. (Ref. 35). 
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Finally the rigorously measured 1200 °C FeCrAlY data of Naumenko et al. (Ref. 35) are presented 
here as open symbols for k P j-G,. Flere grain size data (G) and scale thickness (x) (from Figs. 6 and 10 in 
Ref. 35) were used for samples oxidized for 24, 93, 500, and 2000 h. The quantity dx/df was estimated 
graphically (from curves in Fig. 10 of Ref. 35). Thus k pi = 2x dx/dt and n, = k p j-G I were evaluated from 
the raw experimental data (symbols). It is also noted that IT, may be obtained more analytically from their 
relations for x(t) and G(t) determined from the continuous 100-h isothermal oxidation TGA data. It can be 
shown that 2C = k pj G, = n„ a constant equal to 0.4546 pnrVh ( 1 .263x1 0 4 pnrVs) (dashed line) and 
implies 8D sh for the more general case. 


4.2 FeCrAl(Zr) Oxidation Results (This Study) 

With well-behaved kinetic behavior established under defined microstructural conditions, 
experimental k pj and G, can now be used to predict 8/J„ b derived from scale growth data. Again the 
W-W-S solution (Eq. (10)) is employed, but in reverse. The resulting behavior is given in the Arrhenius 
plot of Figure 9. The same activation energy and correlation coefficient of r 2 = 0.98 are obtained, and the 
full relation given as- 


ST>gb o (oxidation) = 1.79 x 10 10 exp 


-375kJ 


m 3 /s [W-W-S] 


(19) 


y Wada prediction (Eq. (8)), current study 
Regression fit extrapolation 

FeCrAI oxidation: 

0 W-W-S solution (Eq. (10)), current study 

Regression fit extrapolation 

A LNP solution (Eq. (1)), current study 

Regression fit extrapolation 

□ From Naumenko (Ref. 35) 

G} From Messaoudi (Ref. 36) 



Figure 9. — Oxygen grain-boundary diffusion product 8Dgp q estimated 
from FeCrAI(Zr) oxidation and W-W-S solution, Equation (10) of current 
study. Shows agreement with literature value calculations, but offsets 
from permeability predictions. 
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The same offset is seen from the values predicted from permeability relations here as found for k pj in 
Figure 7. Values of 8T) g b.o,mt. determined from the Wada, Matsudaira, and Kitaoka (Ref. 24) relations in 
Equation (1) using the interface values of P 0n eq in Table 1 yield 


i -y 

8D b o (permeability) = 5.16x10 ~ exp 


-298 kJ 
RT 


m 3 /s [Wada, 5/J gh ] 


( 20 ) 


Also, 5Dg b>0 treated as an average or constant value across the P 0i gradient (outside the Wagner 
integral) yields the LNP solution (Eqs. (8) and (11)): 


5Dg bo (° x idation) = 8.12x10 11 


exp 


394 kh 
RT ) 


m 3 /s [LNP] 


( 21 ) 


This is seen in Figure 9 (filled triangles) to be below that produced by using the experimental data in 
the W-W-S solution, again by the factor A lnEg, /6 as expected from Equation (12). The differences 

- 1 / 

reflect the A.P 0 /6 versus A In P () _ integration limits for the two expressions in the Wagner integral. 

Also shown are literature values obtained as in Figure 7, when the grain size was known well enough 
to estimate the oxidation product FI. They are seen to essentially lie on the same curve as the present data. 


4.3 Comparisons to Oxidative Literature Values for 5D* b 0 


Now Figure 10 compares the 5Dg b 0 results for the LNP solution with similar values obtained from 

the oxidation literature. Data from Young et al. (Ref. 12), obtained from the same FeCrAlY alloy as 
Naumenko et al. (Ref. 35), are included as the cross-square symbols (Eq. (16) in Ref. 12). They are seen 
to fall close to values predicted for the present data using this LNP formula (i.e., assuming that 8D gb ,o is 
independent of P Qi ) . Also included are data similarly obtained for YjCE-dispersion-strengthened 

FeCrAl, MA956 (Refs. 13 and 14), with a reported activation energy of 294 kJ/mole, as compared to the 
394 kJ/mole determined from the LNP analyses of the present FeCrAl(Zr) data. The same alloy had been 
characterized earlier by Ramanarayanan et al. (Ref. 10). Using their reported average k p and grain size, 
the LNP results are again seen to fall exactly on the line representing the present study, with an activation 
energy of 388 kJ/mole, and remarkably close to those reported by Clemens et al. for this alloy (Ref. 14). 
However, Ramanarayanan et al. (Ref. 1 0) also performed their own 8D gbi o analyses but arrived at a 
different expression: 


SDg b o (oxidation) = 4.92 x 1 0 8 



— 560 kJ 
RT ) 


m 3 /s 


( 22 ) 


This was seen to diverge to substantially higher values of 8D gbi0 and is not shown here. It is also said 
to be adjusted by P Q 76 , presumably for the environmental test atmosphere, the gas surface. 

Furthermore, an analogous LNP analysis was performed on the MAX compound, Ti XlCa (Ref. 32), 
shown as the crossed diamonds in Figure 10. Except for short times, especially at lower temperatures, the 
scales were primarily a-Al 2 0 3 . Their results are very close to that presented for the four FeCrAl alloys in 
Figure 10, and linear regression produced 375 kJ/mole. In total, the plot in Figure 10 shows reasonable 
agreement from seven separate studies. Although this agreement in experimental data is encouraging, it does 
not necessarily indicate that the LNP solution is the preferred method of calculating 8D gb ,o- It just means 

that there are many studies with an available 8Z)g b 0 determined this way. Had all the studies provided k pj 
and corresponding G, data one could evaluate the same studies via the W-W-S solution as well. 
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— Wada prediction (Eq. (20)), current study 

Regression fit extrapolation 

“ “ W-W-S solution (Eq. (10)), current study 

Regression fit extrapolation 

“ ” LNP solution (Eq. (8)), current study 
Regression fit extrapolation 
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FeCrAI oxidation: 

A Hoskins 875 FeCrAI(Zr), current study 
BB FeCrAlY, Young et al. (Ref. 1 2) 

^ Ti 3 AIC 2 , Wang et al. (Ref. 32) 

MA956 FeCrAI-Y 203 , Clemens et al. (Refs. 13 and 14) 

□ FeCrAlY, from Naumenko et al. (Ref. 38) 

O Imphy FeCrAI(Zr.Si), from Messaoudi, Huntz, and Lesage (Ref. 39) 


MA956 FeCrAI-Y 203 , from Ramanarayanan (Ref. 10) 



Figure 10. — Average oxygen grain-boundary diffusion product 8Dgb q for 
FeCrAI(X) alloys. Literature values from Young (Ref. 12), Wang (Ref. 32), 
and Clemens (Ref. 14) agree with those obtained here by LNP solution 
for alumina scale growth. 
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Regarding the two W-W-S and LNP techniques used to calculate k pj from 5D gb 0 , it is again useful to 
recall the standard factor relating these terms as stated in Equations (8) and (10). Here it was shown that 

the factor relating k pp (LNP), as a constant or average 5D/ b 0 , to the k pj (W-W-S) was A In P () ^ /6 
(Eqs. (11) and (12)). This offset was indicated in Figure 7. Conversely, any &D gb ,o determined from the 
same experimental k p j using a constant 5/J/ b 0 outside the Wagner integral will be the inverse factor of 
6/A In To,, times SY> gbi o,int determined from the W-W-S relation. This factor is so indicated in Figures 9 

and 10. It is thus seen that the LNP analyses predicted k p further from the experimental k p (Fig. 7) or 
SOgb.o.int further from the permeability derived relations than did the W-W-S analyses (Fig. 9). 

4.4 Comparison With 5Z> g b,o Measured by ls O Tracer Studies 

Finally, as with most discussions of diffusion in alumina scales, a more global perspective is provided 
by comparing these predictions with summary plots of diffusion in alumina, such as those recently 
discussed by Heuer et al. (Refs. 22 and 40) in Figure 11. Here 80 gb ,o data measured from 18 0 tracer 
experiments for both bulk alumina (under no P 0i gradient) and alumina scales (growing in a large 

gradient) are compared. The top curve (i) is that obtained by applying the Wada, Matsudaira, and Kitaoka 
(Ref. 24) relation for diffusivity to the low Pq 2 interface on a growing scale (Eqs. (1) and (20), Table II). 

The middle curve (ii) is that predicted from the present data for alumina scales formed on Hoskins 875 
FeCrAl(Zr) (Eqs. (10) and (18)). Here oxidation kinetic and grain size data use the W-W-S solution, 
predicting diffusivity from II, (k pj and G,). The next curve (iii) is that obtained from the same oxidation 

data, but now using the A lnR 0 , solution (LNP, assuming a constant average 8D/ b 0 , Eqs. (8), (11), and 

(21)). The individual data points for curves 1 to 3 had been shown in Figures 9 and 10. The lowest curve 
(iv) was obtained from the permeability-derived equation for SY> gb ,o,gas for the case of exposure to 1 atm 
air, with no P Qn gradient (Eq. (1), Table II, Fig. 3). 

The 18 0 tracer studies for bulk alumina refer to pure or Mg-doped AFO3 for Reddy and pure or Y-doped 
AI2O3 from Prot et al. (Refs. 41 to 43). Only curve 4, appropriate for bulk alumina without a P 0i gradient, 

appears to approach magnitudes of the undoped or Mg-doped AFO3, but the data are quite high compared to 
that extrapolated from Y-doped A1 2 0 3 . Curve 4 does not have an activation energy similar to any of these, 
but at least is much closer than any predictions for SO gb ,o appropriate for oxidation (curves 1 to 3). 

The ls O tracer studies for double oxidation of FeCrAl alloys refer to those of Messaoudi et al. 

(circles, squares, and diamonds) (Ref. 36) and Chevalier et al. (triangles) (Ref. 44). The 1 100 °C point for 
MA956 (Refs. 14 and 45) is also included (dotted hexagon) as it corresponds to fitting an 18 0 SIMS 
profile. All the alloys of Messaoudi are in the vicinity of curve (iv) for bulk alumina, whereas the 
undoped alloy of Chevalier (open triangles) is grouped near curve (ii) predicted for growing scales. The 
Y-doped alloy (filled triangle), however, is near the Messaoudi data and curve (iv). The Clemens/ 
Bongartz point falls directly on the permeability LNP solution, as in Figure 10. Thus only some 18 0 tracer 
data seems comparable to the appropriate permeability predictions. These inconsistencies are difficult to 
resolve and underscore the challenges encountered with measuring or predicting diffusivity in alumina 
scales and comparing the same to that measured for bulk alumina. (It should be noted that more 
measurements for bulk alumina from Reference 44 (not shown) fall between the oxidation values, far 
from other bulk data, and are not reconciled here). Overall it is seen, as in Figure 3, that the P 0i gradient 

in a scale can result in diffusivities expected to be 3 to 5 orders of magnitude greater than for bulk 
alumina in no P Qi gradient. It also appears that some double oxidation experiments may reflect this in 

varying degrees. 
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Predicted: 

— Wada, for growing scales 
Regression fit extrapolation 

— W-W-S solution, current study 

— 8Dgb o LNP solution, current study 

— Wada, for bulk (in air, no Pq 2 gradient) 

Regression fit extrapolation 

18 Q Tracer Studies: 

O Bulk Al 2 0 3 , Reddy (Ref. 41 ) 

© Bulk doped with Mg, Reddy (Ref. 41 ) 

Regression fit 

Bulk, Prot (Ref. 42) 

Bulk with Y, Le Gall (Ref. 43) 

© FeCrAI(Zr) scales, Messauodi (Ref. 39) 
ffl FeCrAI(Zr.Y) scales, Messauodi (Ref. 39) 

© FeCrAI(Y) scales, Messauodi (Ref. 39) 

Regression fit 

A FeCrAI scales, Chevalier (Ref. 44) 

A FeCrAI(Y) scales, Chevalier (Ref. 44) 

Regression fit 

© FeCrAI scales, Clemens (Ref. 14) and Bongartz (Ref. 45) 



Oxidation temperature, 1/T, 10 4 K 1 

Figure 1 1 . — Comparison of oxygen grain-boundary diffusion product 5D g b t o 
calculated from four permeability relations (solid lines) to those measured 
from 180 tracer profiles (dashed lines) in bulk Al 2 0 3 and AI 2 Q 3 scales. 
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Finally, comparison is now made to more recent data for individual grain boundaries in bulk alumina 
(Fig. 12). The lower band (oriented) refers to bicrystal studies with varying habit planes and coincident 
site lattice I values (Ref. 46). The values for two 17 and S3 1 boundaries are quite close to those projected 
from the permeability relations for bulk alumina at atmospheric conditions, whereas the remaining curves 
are well below these. The upper cluster of points (random), recently obtained for individual grain 
boundaries in bulk polycrystalline (Ref. 47, personal communication), appears in a band with values near 
the projection of the W-W-S curve for scales. Since these represent boundaries where substantial tracer 
concentrations were observed, they may represent an upper band of diffusivities. In all, a variability of up 
to 8 orders of magnitude may be concluded from the bulk alumina data. Clearly, the nature of each grain 
boundary is of great importance and will enter into any discussion of alumina scale growth. At present, 
we can only assume that scale growth is controlled by some distribution of countless grain-boundary 
orientations, perhaps affected by growth textures in special cases. 

Predicted: 

— Wada, for growing scales 

Regression fit extrapolation 

H W-W-S solution, current study 
“ — - SDgb.O LNP solution, current study 
— Wada, for bulk (in air, no Pq 2 gradient) 

Regression fit extrapolation 

Measured: 

O Individual random boundaries, Nakagawa and Heuer (Ref. 47) 

Oriented boundaries, Nakagawa et al. (Ref. 46): 

Line X hkil 

a 31 71140 

b 7 2310 

c 21 2310 

d 21 4510 

e 7 4510 



Figure 12. — Comparison of oxygen grain-boundary diffusion product 5D g b,o 
for bulk Al 2 0 3 and Al 2 0 3 scales calculated from the four permeability relations 
(see Fig. 11) and those measured by 18 0 tracer studies of individual grain 
boundaries: (1) Random boundaries measured by ToF-SIMS (time-of-flight 
secondary ion mass spectrometry) profiles of polycrystalline alumina (Nakagawa 
and Heuer, Ref. 47, personal communication) and (2) boundaries from bicrystals 
oriented by habit plane and tilt angle (Nakagawa et at, Ref. 46). 
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5.0 Activation Energies 

A catalogue of activation energies obtained is presented in Table V. The top entry indicates the 
thermodynamic parameter resulting from an Arrhenius relation of the equilibrium P 0 ^ for the oxidation 

reaction (Fig. 2). The temperature dependence of the free energy of AI 2 O 3 formation provides the 
dominant effect here. The next two values are those determined by Wada, Matsudaira, and Kitaoka 
(Ref. 24) for oxygen and aluminum boundary diffusion under atmospheric conditions, with P 0i 

~0.2 atm, independent of temperature (Fig. 2). The following two apply to the scale-metal interface 
(Table 11 and Fig. 2). It is seen that Q mt = (9 perm + nQ e q , where (Q ml is the activation energy resulting for 
diffusivity at the P Qi calculated for the scale-metal interface, Q pcm] is that provided from the 

permeability-derived diffusion relations (Eq. (1)), n is the exponent of the Pq term in the latter, and Q eq 

is that produced by the temperature dependence of P 0i eq . Thus, the activation energies for oxygen and 

aluminum diffusivity at the interface are reduced by 1/6 x (1012 kJ/mole) and increased by 3/16 x 
(1012 kJ/mole), respectively. 

TABLE V.— ACTIVATION ENERGIES REPORTED FOR oD gbo AND 6D gbAi 


[Calculated from permeability equations and estimated 8Z) g b,o for FeCrAl(X) alloys using 
Wagner- Wada-Smialek (W-W-S) and AlnP 02 solutions.] 


Source 

Application 

Value 3 

Activation energy, 

Q . 

kJ/mole 

Reference 

Thennodynamic 

A1 2 0 3 - FeCrAl 

P 

r 0 2 ,eq 

1012 

Jacobson (current study) 

Wada 

Bulk, scale (gas) 

SDgb.o 

467 

Wada (Ref. 24) 

Wada 

Bulk, scale (gas) 

SDgb.Ai 

604 

Wada (Ref. 24) 

Wagner- Wada 

Scale, interface 

Sf^gb.O 

298 

Current study 

Wagner- Wada 

Scale, interface 

8^gb,Al 

794 

Current study 

FeCrAl(Zr) 

Scale (measured) 

kp,i 

375 

Current study 

FeCrAl(Zr) 

Scale (W-W-S) 

8^2gb,o 

375 

Current study 

FeCrAl(Zr) 

Scale (LNP) 

8^2gb,o 

394 

Current study 

FeCrAl-Y 2 0 3 

Scale (W-W-S) 

8^2gb,o 

388 

Ramanarayanan (Ref. 10) 

FeCrAl-Y 2 0 3 

Scale (LNP) 

V)gb,0 

294 

Clemens (Refs. 13 and 14) 

FeCrAlY 

Scale (LNP) 

SDgb.o 

321 

Young (Ref. 12) 

Ti 3 AlC 2 

Scale (LNP) 

SDgb.o 

375 

Wang (Ref. 32) 

NiAl(Zr) 

Scale (creep) 

D g b 

434 

Veal (Ref. 48) 


d ^c> 2 , eq is equilibrium oxygen pressure, 8D g b is grain-boundary diffusion product (for oxygen, O, and 
aluminum, Al), and D g b is grain-boundary diffusivity. 


The k pq experimental data of the present study are given next and yield 375 kJ/mole. This is 
maintained in the Wagner-Wada-Smialek derivation of SjDg b ,o, but changes to 394 kJ/mole in the A In Pq 

analysis, similar to 388 kJ/mole determined from the Ramanarayanan data. The last three list A In P () _ 

results, giving 294, 321, and 375 kJ/mole for 5Z ) gbj0 from oxidation of MA956, FeCrAlY, and Ti 3 AlC 2 , 
respectively. The first value (Ref. 14) may have been influenced by the high value at 900 °C; the next 
series (Ref. 12) was obtained from a narrow 100 °C temperature range, and the last (Ref. 32) is in closer 
agreement with the present study. All of these are substantially lower than the bulk alumina studies shown 
in Figure 1 1 . 
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6.0 Additional Transport Issues 

Up to this point the coupling of the permeability relations with oxidation kinetics has concentrated on 
oxygen grain-boundary diffusion in pure alumina. To some extent, this is all that is practical given that it 
is difficult to reliably separate any A1 outward growth from primarily inward O growth. Furthermore, the 
study of dopant effects on diffusivity in bulk alumina is ongoing, with perhaps a tenuous connection to 
thermally grown scales with various grain-boundary dopant distributions. 


6.1 Aluminum Diffusion 


The importance of aluminum outward diffusion is often discussed in the literature, but it is not 
considered here given the detail required to cover just oxygen diffusion. For the sake of completeness, the 
same Wagner-type integrations were performed for the Wada relation for 8Z3 g b,Ai (Eq.(2)) and are listed in 
Appendix C. It is noteworthy that similar reduced terms result, allowing the A1 diffusional component of 
k p to be given as a simple function of SD gb ,Ai,gas (Eq. (C4)), such that 


k 


p,i, A1 


16 8D 


W-W-S 

gb.Al.gas 


G, 


(C4) 


and the W-W-S expression is related to the LNP solution or average diffusivity by the factor 3/16 
A In Eg, . The difficulty, as with any system exhibiting mixed transport, is experimentally differentiating 

the minor A1 contribution from that of oxygen to k pj (Refs. 29 to 3 1). 

We expect that the majority of transport is by oxygen diffusion, as controlled by the interface P Qi , 

such that oxygen diffusivity at the interface is 2 to 5 orders of magnitude greater than aluminum at the gas 
surface, the largest differential being at lower temperature. Furthermore in the modeled case of oxidation 
at 1100 °C, A1 diffusivity is seen to be negligible except at the outermost near surface of the scale (Fig. 10 
in Ref. 18). Thus the focus on oxygen diffusivity in the present study is consistent with the permeability 
projections. 


6.2 Dopant Effects 

It must be acknowledged that A1 transport is a key factor in explaining growth stress and the 
formation of oxide within the scale or at outer grain-boundary ridges. The effects of impurities or 
intentional dopants often enter this conversation, as in producing charge compensating defects that may 
affect A1 diffusion. Although, again, not the focus of the present study, the permeability experiments 
performed on Lu-, Y-, and Hf-doped alumina addressed these issues, reaching P 0i down to 10 s Pa 

(Refs. 17, 18, and 49). It was found in general that Lu and Y dopants (on the order of 0.4 at.% LU 2 O 3 ) 
reduced oxygen grain-boundary diffusivity by about 3 times, with no change in aluminum diffusivity, 
whereas Hf dopants (0.2% Hf 02 ) reduced A1 diffusivity by about 2 times, with no change in oxygen 
diffusivity. This is supported by an order of magnitude reduction in oxygen diffusivity for an Y -doped 
bicrystal found in tracer studies (Ref. 50). Surprisingly, co-doping with equal amounts of Lu and Y 
produced about a 1.5-times increase in both A1 and O diffusivity in the permeability studies (Ref. 18). 
These results support the often-observed beneficial effect of Hf dopants on alumina scale growth, but do 
not completely resolve effects typically reported for optimum Y-doping or co-doping where A1 diffusivity 
and outward growth appear to be reduced. In the present study, it may be expected that some of the 1 0 to 
40 times lower offset of oxidation-derived values in comparison to permeability predictions (Figs. 7, 9, 
and 10) are due to dopant effects. In this work and the Messaoudi data (Ref. 36), the dopant was Zr and 
may be expected to exhibit effects similar to those of Hf. However, examination of these figures shows 
that the data for Zr, Y, Ti, and Y 2 0 3 dopants all agree with nearly a single Arrhenius relationship. Thus 
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clear discrimination of oxidation kinetics by dopant type, to the extent suggested by the permeability 
studies, has not emerged. 

There has been much discussion of various effects of Y-doping (or other reactive elements). Initially 
reduced A1 transport, growth stress, and buckling occur as Y levels increase and then give way to 
increased weight gains and disruptions from oxidation of Y-rich precipitates (e.g., YAG, yttrium 
aluminum garnet) and oxide intrusions. Finally, Hf dopants are recognized as perhaps the most effective 
in reducing k p , by about an order of magnitude, and receive a great deal of attention in protective 
commercial Hf-modified MCrAlY overlay or NiAl diffusion coatings (Ref. 4). There does not yet seem to 
be an exact correlation between dopant effects on oxidation with those predicted by the permeability 
studies. 

It is acknowledged that 10 to 100 ppm of common Ni, Co, Fe, Cr, and Si base metal impurities may 
be present in alumina scales grown on typical alloys. A comparison of reactive element dopant effects in 
bulk aluminas containing these metals may thus provide further insight if data from a “pure” alumina 
scale were analyzed, such as that formed on Pt-Al alloys. 

6.3 IT Decay for FeCrAl(Zr) 

According to Equation (10) the oxidation product n = k pj - G, should be a constant and equal to 12 
SDgb.o This was indeed the case for the 1200 °C FeCrAlY data in Figure 8 (Ref. 35), even out to 2000 h. 
However, the present study suggests an accelerating decrease in this term after about 200 h. One potential 
cause might be reduced grain-boundary diffusion according to Equation (10). Although there is no 
evidence either way, significant A1 depletion would reduce the interfacial A1 activity and increase P (u eq 

according to Equation (14). Then an increase in P (u would result in a decrease in 8E> gb , 0 according to 

Equation (1). Since the A1 content and thickness of the Naumenko samples were less than the present 
FeCrAl(Zr) alloy, it cannot be argued that greater A1 depletion from a reduced A1 reservoir resulted in the 
present study. Nor was greater spallation observed. At present there is no compelling explanation for the 
long-term decay in the oxidation product, IT. 


7.0 Summary and Concluding Remarks 

The temperature and pressure dependence of oxygen and aluminum grain-boundary diffusivity were 
determined by permeability studies of pure bulk alumina wafers by Wada, Matsudaira, and Kitaoka in 
2011. Their work is extended in the present study to further elucidate the growth kinetics of high- 
temperature alumina scales as a function of grain-boundary diffusivity, 8D gb . First, equilibrium 
thermodynamics was applied to the case of a commercial FeCrAl alloy to determine the equilibrium 
oxygen pressure at the metal interface P 0i int eq . These values and the permeability relations were 

applied to a modified Wagner equation for ionic scales. For growth primarily by oxygen grain-boundary 
diffusion, the following relation was developed (Eq. (18)) for the instantaneous parabolic rate constant 
(at a grain size of 0.5 pm): 


k i [W- W - S] - 1 .24 x 1 0“ 4 expf 298 ^ 

V RT 


m 2 /s 


This is compared to a relation determined experimentally (exptl) for the same scales (Eq. (16)), with a 
moderately higher activation energy: 


i', exptl = 4.29 x !0 3 exp^- 


-375kJ 

RT 


m 2 /s 
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It is also found that a grain-size invariant oxidation rate product, II = k p y G t is a simple function of 
oxygen grain-boundary diffusivity from a Wagner-Wada-Smialek (W-W-S) derivation. Thus, from 
Equation (10), 


- k pi Gj - 1 25/J gbX) int 

here 5/J gb is a strong function of P 0 , and varies across the scale. This provides a surprisingly convenient 

means of drawing comparisons: it allows either prediction of oxidation rates from independently measured 
diffusivity or prediction of grain-boundary diffusivity from instantaneous oxidation rate and grain size. As 
simple as this relation appears, the real complexity of transport issues are embedded in the starting 
permeability relations, the Wagner integral, and the Hart/Smeltzer diffusivity equations. 

Comparisons to similar oxidation studies yielded relative overall agreement, confirming some level of 
confidence in the experimental data and the concept of an invariant oxidation constant, II. The predicted 
oxidation kinetics trended from ~10 to 40 times greater than the experimental values as temperature was 
reduced. Although this is not exact agreement, it is enormously closer than extrapolations of bulk alumina 
diffusivity studies obtained with no P 0i gradient. Moreover, the discrepancy may be partially ascribed to 

the experimental data pertaining to doped alumina scales, whereas the predicted data are for pure alumina. 
For the sake of completeness, comparisons were made to predictions where constant (average) 

diffusivity was assumed across the scale. Here a standard correction factor, (a In P ()n )j6 could be 

applied to obtain k pj from the comparable W-W-S expression or the inverse factor for obtaining 
diffusivity from the experimental k p j. Also, the extremes of aluminum and oxygen diffusivity were 
explored at each interface, and SD gb ,Ai was found to be negligible compared to ST> gb ,o, except at the gas 
surface. Though not utilized here, an analogous method resulted in a simplified W-W-S relation for 
8D gh ai, which can just be added to the oxygen contribution for a full solution to the Wagner integral 
(Eq.(B4)): 


k p ,i,MGi = 168 D 


W-W-S 

gb,Al,gas 


In all, the permeability studies, by quantifying P Qi effects on grain-boundary diffusivity, have allowed 

for improved insights and closer connections between predicted diffusivity and alumina scale growth 
kinetics. 
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Appendix A. — Symbols 


Symbols: 

A pre-exponential exponent from Equation (9) 
a activity 

D diffusivity 

/ area fraction of short-circuit grain-boundary paths 
G grain size (diameter) 

k oxidation rate constant 

K equilibrium constant 

P pressure 

Q activation energy 

R ideal gas constant 

r 2 correlation coefficient 

T temperature 

t time 

x scale weight gain, or equivalent thickness 

z ionic charge 

8 grain-boundary width 

If oxidation product 

Subscripts: 

0.5 0.5-pm grain size 

I low- p o 2 surface 

II high- P ()i surface 

a anion 

A1 aluminum 

c cation 

eff effective 

eq equilibrium 

exptl experimental 

gas gas 

gb grain boundary 

1 instantaneous 

int interface 

L lattice 

m grain growth time exponent 

n oxidation time exponent 

O oxygen 

0 2 oxygen gas 

p parabolic 

perm permeability 
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Appendix B. — Defect Equilibria, Electroneutrality, and P Qi 

The oxidation reaction and basic defect equilibria for undoped alumina can be described by Equations 
(Bl) to (B5) using Kroger-Vink notation (Ref. 23). Equations (B2) and (B4) apply to the oxygen vacancy 
concentration [Fq*J at the low- P () surface (n-type), and Equations (B3) and (B5) apply to the aluminum 
vacancy concentration [V"\ ] at the high- P () surface (p-type). The equilibrium constants are then written 

as Equations (B6) to (B8), corresponding to Equations (Bl) to (B3). Solving the electroneutrality 
relations Equations (B4) and (B5) for [V,] and substituting into the corresponding equilibrium relations 
Equations (B7) and (B8) produce Equations (B9) and (B 11), respectively. Finally, combining Equation 

(Bl) with (B9) produces Equation (BIO). Thus, at the low- P () surface [f o ’‘] is proportional to P c ''j , and 
at the high P () surface [Vq] is proportional to . 


2Al + %0 2 = A1 2 0 3 

(Bl) 

2A1 = 2 Alf j + 3 Vq + 6 e' (n - type) 

(B2) 

%0 2 =30§ + 2TXi +6/?* (p-type) 

(B3) 

4'o]=[e'} 

(B4) 

3[^i ]=k 

(B5) 

Kl= 

(B6) 

i 2 klVr 

^2 — L 2 

a A1 

(B7) 

,, iwkf 

A 3 - 3 / 

p/2 

o 2 

(B8) 

is _ z v O 
K 2~ 2 

a Al 

(B9) 

Vq =2~ 2 A K A P~ l//& 

0 U.J ° ! 

(BIO) 

(74= 3-3/4 TAf 16 

(Bll) 
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Appendix C. — Relations Derived for Aluminum 
Grain-Boundary Diffusivity, 8Z) gb?A 1 

The complete Wagner expression for alumina growth by aluminum and oxygen diffusion is obtained 
from Equation (4) as 


'02,gas 


k„ ; = 2x— - f (3 /d 


'p,l 


d t 


+ D, 


2 ^eff, A1 ^ ^eff,0 


■f,o) 


din P ( 


O, 


M3 2 ,int 


(Cl) 


Dealing again only with one component, the aluminum contribution, obtained from Equation (2) 
(repeated here): 




= 2.475 xio -5 



-604kJ V, + & 
RT ) °2’ n 


produces the following form, Equation (C2), integrated to Equation (C3) (Wagner-Wada-Smialek, 
W-W-S) and simplified to Equation (C4) (and (C5) if a constant Z)* bA1 is assumed). Finally, the 
relationship between the W-W-S and LNP solutions for 5/J gh ,Ai are given by Equation (C6): 


f p,i, A1 


mC> 2 ,gas 


r 02 ,int 


2(5/4 


gb,Al 


G; 


dlnP 


kp,i, A1 


<2 ex P( ~ 2 R T 


1 6AP+^ 6 
u 2 


1 65D 


f p,i, A1 - ' 


W-W-S 

gb,Al,gas 


G, 


(C2) 


(C3) 


(C4) 


k 


p,i, A1 


C ; 


AlnPo 2 


5Z> 


W-W-S 

gb,Al,gas 


3A In P 0i 
16 


S^gb.Al 


(C5) 

(C6) 
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